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20  abstract  ''Clinffnue  on  r«i  fr$a  ntda  If  neca<itary  and  tdantifv  blork  nnnyh^r 

.'I  he  field  strength  of  the  l■om^ler(■ial  radio  hroadcas;  station  KHI.E  was 
measured  at  sevi'ral  hundred  locations  in  the  desert  .soullieasi  of  ^■uma, 
•Ar'zona.to  study  ME  ground-wave  propagalion  oxer  tin-  desert  \  allevs  and 
mountain  ridges.  Field  strengths  below  the  surfac'  of  the  desert  iloor  were 
measured  in  .speciallv  drilled  holes  to  observe'  the  penetration  of  M  F  into  the 
desert  soil.  Idealized  propagation  patn  models  are  shown  to  produce  effects 
similar  to  the  observations. 
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MF  Radio  Field  Strength  Measurements 
in  Desert  Terrain  Near  Yuma,  Arizona 


1.  IM'RODUCTION 

The  possibility  that  a  network  of  underground  tunnels  or  silos  containing  the 
MX  mobile  missile  might  be  constructed  in  the  desert  regions  of  the  southwestern 
United  States  has  renewed  interest  in  techniques  for  survivable  communications 
between  buried  terminals  which  may  be  up  to  50  km  apart  in  adjacent  valleys  sepa¬ 
rated  bv  high  mountain  ridges.  One  proposed  method  is  to  use  an  MF  radio  link 
with  the  transmitting  and  receiving  antennas  located  on  the  underground  mobile 
launch  control  centers.  The  propagation  path  for  such  a  link  would  include  two  seg¬ 
ments  through  the  trench  wall  and  soil  overburden  in  addition  to  that  over  the  desert 

terrain.  Although  MF  propagation  curves  over  homogeneous  earth  are  readily 
1  -4 

available,  propagation  anomalies  caused  by  topographical  and  electrical  irregu¬ 
larities  unique  to  the  proposed  MX  sites  need  to  be  studied  in  situ. 


(Received  for  publication  24  November  1980) 
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2.  EXPERIMENTAL  PROCEDURE 

Between  11  September  and  24  September  1977  personnel  from  Rome  Air 
Development  Center  (RADC/FEP)  made  field  strength  measi’rements  of  the  Yuma, 
Arizona  broadcast  station  KBLU  (560  kHz,  1  kW,  omni -di r  ectional  in  the  daytime) 
at  several  hundred  positions  on  Luke  A  I’  Bombing  and  Gunnery  Range,  a  desert 
region  southeast  of  Yuma.  At  each  site  the  surface  magnetic  field  was  measured 
using  a  Radio  Interference-Field  Intensity  Measuring  Equipment,  Singer  Model 
NM-25T,  with  a  shielded  15-in.  -diameter  loop  antenna  oriented  for  maximum 
signal.  At  some  of  the  sites  the  loop  also  was  lowered  2  m  into  a  hole  drilled  by 
a  power  auger  to  probe  the  radio  wave  attenuation  with  depth.  The  measurements 
were  conducted  only  during  daylight  hours  to  eliminate  the  necessity  of  correcting 
for  KBLU  day-to-night  antenna  pattern  changes  and  to  minimize  sky  wave  inter¬ 
ference. 

The  measurement  sites  were  accurately  located  in  almost  every  case  bv 
choosing  positions  near  the  numerous  USGS  markers  found  alongside  Jeep  trails. 

The  sites  are  shown  in  Figure  1,  along  with  the  location  of  the  KBLU  transmitting 
antenna  in  downtown  Yuma.  The  distribution  of  the  markers  was  such  that  no  series 
of  measurements  was  strictly  aligned  along  a  radial  path  from  KBLU.  However, 
six  sectors  were  chosen  within  which  the  terrain  could  be  characterized  reasonably 
well  by  a  single  profile,  thus  providing  a  convenient  wav  to  group  and  display  the 
data.  Dashed  linos  are  used  in  Figure  1  to  show  the  sector  boundtiries.  A  repre¬ 
sentative  path  profile,  corrected  for  the  stundtird  tropospheric  refractive  index 
gradient^  (4/2  earth  radius),  was  constructed  for  each  sector  using  USGS  maps. 

3.  FLAT  EARTH  PROPAGATION  THEORY 

Over  a  flat,  homogeneous,  and  well-conducting  earth  the  vertical  electric  field 
St rength-dista.nce  product  is  giver,  bv 

■E-xi  =  300  }/p-  !F(p^)|  (Volts)  (1) 

where  F  is  the  effective  radiated  pow'er  in  kW, 

F(Pp)  =  1  -  i  ^77  p^  *0  ^  •  erfc(i  ^ p^ ) 

a-  Reference  Utita  for  Radio  Engineers,  5th  ed.  ,  ( 1957)  ITT  Corp.  ,  New  A  ork 
p  74  1  ff.  ° 


and 
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F(p  )  is  the  Somnierfeld  propagation  factor,  p  is  the  numerical  distance, 

®  - 1 2  ® 
t  =  8.  854  X  10  F/m,  u;  is  the  radian  frequency,  CT  is  the  surface  conductivity 

(Siemens /meter)  and  kx  is  the  distance  in  wavelengths.  Eq.  (1)  can  also  be  expressed 

log  Ie-xI  =  log  IfI  +  log  300  yp  .  (2) 

where  log  =  logjg.  1  f|  is  plotted  us  a  function  of  p^  on  a  logarithmic  scale  in 
Figure  2.  For  numerical  distances  of  the  order  of  unity  or  less, 

log  I  fS  :  -0.  183  tppl  .  (3) 

At  560  kHz,  Fq.  (3)  is  valid  for  distances  out  to  at  least  50  km  if  ct  =  0.  01  S/m  and 
to  100  km  foro  =  0.02  S/m.  Combining  Eqs.  (2)  and  (3), 

0.  Of)  13  £  a;k 

log  'e  -  x!  -  -  - g— 2 -  X  f  log  300  yp  .  (4) 

Equation  (4)  is  in  the  standard  form  y  =  mx  +  b,  where,  upon  the  application  of 
linear  regrcs.sion  analysis  to  the  measured  (E'  x(  data,  the  slope 

t  ol'  k 

m  =  -0.  0915  —  (5) 


yield.san  estimate'  (or  u,  and  the  y-intercept 

b  =  log  kOoVi*  (6) 

indicates  an  apparent  radiated  power 


F'igure  2.  Amplitude  of  the 
Propagation  Factor 
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Product 


I'igure  6.  Coyote  Pea'K 
(110“’  -  116  =  )  Sector  Path 
Profile  (lower  portion)  and 
Signal  Strength-Distance 
Product 
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4.  RKSII LTS  OF  REGRESSION  ANALYSIS 


4.1  Estimate  uf  Conductivity 

The  upper  portions  of  Figures  3  to  8  show  measured  values  of  |  F  •  xl  on  a 
logarithmic  scale  vs  distance  from  the  KBLU  antenna.  The  lower  portions  are 
corresponding  terrain  profiles  which  reveal  long  flat  path  segments  occasionally 
interrupted  by  relatively  steep  ridges.  The  effective  conductivity  of  these  flat  seg¬ 
ments  was  estimated  by  fitting  a  regression  line  to  selected  data  points  not  too  close 
to  the  mountain  ridges,  and  then  using  Fq.  (5)  to  estimate  o.  The  regression  line 
is  shown  in  the  figures  as  a  solid  line,  and  the  o  estimates  are  given  in  Table  1. 
Assuming  the  data  is  normally  distributed  about  the  regression  line,  we  cun  assign 
confidence  limits  to  the  calculated  values  of  u.  The  results  of  choosing  a  95  percent 
interval  are  shown  in  the  last  column  of  Table  1.  Fvidently  only  the  100°  -  107.  5° 
sector  contained  enough  data  over  a  sufficiently  long  and  flat  range  to  permit  a 
reasonably  accurate  estimate  for  u,  which  was  0.022  S/m.  This  value  is  somewhat 
higher  than  the  0.  01  S/m  estimated  by  Morgan'  or  the  0.  008  S/m  estimated  bv  the 
FCC.  ^ 

Some  of  the  data  scatter  within  the  flat  segments  was  due  to  the  fact  that  the 
measurements  were  not  along  radial  paths,  and  each  data  point  within  a  given 
sector  actually  represents  propagation  over  a  slightly  different  path.  .Additional 
perturbations  were  produced  by  local  anomalies  such  as  ore  deposits,  faults,  power 
and  telephone  lines,  and  railroads,  as  well  as  by  focussing,  diffraction,  and  multi¬ 
path  effects. 


Table  1.  F.stimaled  Fffeclive  Conductivity 


Sector 
[ Degrees  | 

No.  Data 
Points 

.. 

Correlation 

Coefficient 

Fst  CT 
[S/m  1 

05%  Conf.  Limits 
[S/m ) 

90  -  100 

24 

0.  39 

0.  023 

* 

100  -  107.5 

44 

0.70 

0.  022 

0.  016  -  0.  033 

107. 5  -  no 

19 

0.  42 

0.  029 

* 

110  -  116 

21 

0.22 

0.  026 

116  -  120 

6 

0.  22 

O 

O 

o 

123  -  137 

20 

0.  52 

0.  023 

0.  on  -  0.  122 

The  data  scatter  was  too  great  to  permit  an  accurate  determination  of  a  . 


(Uue  to  the  number  of  references  cited  .above,  they  will  not  be  listed  here.  See 
Ri'lerences,  page  20.  ) 
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■1.2  .Apparent  Kadiaterl  I’liwer 

in  I'igur'c  s  A  to  H  the  du.shed  lines  plot  I  1:1  •  xj  from  I'lq.  (4)  for  a  Hat  i-.irth  of 
eonduetivit.v  0.02  S/m  from  a  aiiO  kllz  transmitter  i-adiatirig  1  kU  .  I{(■gl'l  s.sion 
analvsis  on  K.AlfC'  me;i.sured  data  over  path  .segments  where  no  mountain  i-irlges  .ii  e 
interpo.sed  between  the  transmitter  and  observer  (for  exampk’,  1-  igure.s  7  and  h) 
yields  y -inte reept .s  di.splaced  appri',.ximately  l.h  dli  below  tlie  equivalent  1  k\\  nomie.al 
\aliie.  Thi.s  di.serepancy  could  be  due  to  the  \'jM-2a  r  calibration,  sinci-  J.  ■  d..t,! 

supplied  by  KBLL:  is  consistently  higher  than  the  HA!)('  measurement.s.  (in  thcothei- 
hand,  beyond  the  Gila  Alountains  (l-'igures  3  through  6)  the  regression  lines  are  dis¬ 
placed  from  3  to  .a  dll  below  the  nominal  1  k\V  value.  This  reduction  in  apparent 
radiated  power  can  be  interpieted  as  the  effect  of  the  mountain  ridgi'  in  extracting 
energy  from  the  propagating  ground  wave. 

The  losses  may  be  due  to  scattering  by  the  terrain  elevation  changes  and/or 
to  increased  absorption  bv  a  poorly  conducting  path  segment,  but  neither  mechanism 
can  be  adequately  treated  by  flat  earth  theory.  lnde<-d  for  certain  path.s  the  .signal 
stri'ngth  drops  by  as  much  as  15  dli  behind  the  mount.ain  ridges  (1-  igures  (>  and  i)  .md 
then  exhibits  a  recovery  effect  with  distance.  Such  a  variation  requires  the  more 
compreliensivo  theory  which  follow's  in  the  next  section. 


.->.  HROPAG  ATIO.N  OVER  ROUGH.  NON-IIOMOGEAEOUS 1  ERR  AIN 
.5.1  Hiil't'or<r<i  Intogrul  Eipialion 

Calculations  of  the  effects  of  irregular  terrain  can  be  made  via  a  one-dirnen- 

0 

sional  integral  equation  derived  by  Ilufford.  '  His  liq.  (II)  in  the  notation  of 
ITgurc  0  is 


W(x)  - 


f  W(s)-  f(s)  - 
o  ^s(x-s) 


(7) 


where  \S  is  Ihe  atti'nuation  function  for  iht'  Ib-rlz  iiotmiiial.  In  terms  of  the  vertical 
electric  field,  an  a[)f)roximale  cx()rcssion  for  \V  is 


0.  Hufford,  G.  A.  (l')52)  .An  integral  equation  approacli  lo  the  problem  of  wave 
propagation  over  an  irregular  u-rrain,  Quart.  .  .A(ipl.  Math.  it:3!t  1-404. 


1.5 


provided  W  is  slowly  varying.  'I'he  quantity 


l'(s)  =  e 


-  i  tr  /4 


kT  r,  ^'■2] 

2T  *  sttJ 


where  5  is  the  (normalized)  surfaee  impedance 
-in  /  4 


I- 


Figure  9.  Geometry  lor  the  Integral  liquation, 
s  and  X  are  the  horizontal  distances  of  the 
scattering  point  Q  and  the  observer  P, 
respectively,  from  the  source  T.  (v  is.  the  slope 
of  r2,  p  i.s  the  terrain  slope  at  Q,  and  n  is  a 
unit  vector  at  Q  normal  to  the  terrain 


The  factor  takes  into  account  the  terrain  slope  and  the  aspect  angle  with 

respect  to  the  observer.  As  may  be  seen  from  Figure  9, 

=  sin  (p  -  (v) 


where  /3  is  the  terrain  slope  and  a  is  the  slope  of  r,^. 

F'or  gently  undulating  terrain  with  gradual  inhomogeneitio.s,  the  functions  W(s) 
and  f(s)  are  slowly  varying.  To  the  extent  that  VV  and  f  <'an  be  represented  as  con¬ 
stants  within  arbitrarily  small  intervals  s  =  s  -  s  ..an  approximate  solution 

n  n  "  1 

ot  Kq.  (7 )  is 


W 


n  + 


W  f  I 
n  n  n 


(8) 


where  W  is  a  constant  value  for  W’  over  the  nth  interval,  f  is  a  constant  value 
n  n 

for  f  over  the  same  interval  and 
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5.2  Numerical  Solution  for  Special  Cases 

Equation  (8)  was  solved  numerically  via  a  CUC  6600  computer  for  several 
case.s.  In  Figure  10  the  results  of  including  a  Gaussian-shaped  ridge 

GU)  =  h  e.'cp  -9  )^j  (9) 

on  an  otherwise  flat,  homogeneous  earth  are  shown  for  I  vo  values  of  conductivity. 

In  Eq.  (9),  G(x)  is  the  terrain  elevation,  h  is  the  ridge  height  ai  x  =  b  and  w  is  the 
ridge  width  measured  at  G(x)  -  h/10.  This  form  corresponds  to  the  model  chosen 
by  Berrv.  The  attenuation  function  lw|  at  first  decreases  at  the  flat  earth  rate, 
and  then  increases  to  a  ma.ximurn  just  before  the  crest  of  the  ridge.  Behind  the 
ridge  there  is  a  minimum,  and  then  a  partial  recovery  with  increasing  distance 
from  the  source. 

Figure  11  illustrates  the  perturbation  in  lw|  caused  by  a  6  km  wide  segment 
of  conductivitv  ff^,  centered  at  30  km,  in  an  otherwise  homogeneous  flat  path  of 
conductivitv  0.  02  S/m.  Such  a  path  represents  the  much  poorer  conductivity  ex¬ 
pected  in  the  mountain  ridges.  lw|  follows  the  flat  earth  attenuation  rate  out  to 
the  beginning  of  the  segment,  and  then,  depending  on  the  selected  value  for 
decreases  more  or  less  rapidly  until  the  far  edge  of  the  segment  is  reached.  As  in 
the  case  of  the  ridge,  |  W  |  partially  recovers  as  the  distance  from  the  inhomo¬ 
geneity  increases.  (The  abrupt  conductivity  change  in  this  model  does  not  satisfy 
the  conditions  under  which  Eq.  (8)  was  derived,  so  the  results  are  not  valid  inthe 
immediate  vicinity  of  the  edges  of  the  segment.  )  A  comparison  of  Figures  10  and  11 
with  7  and  8  seems  to  show  that  the  effect  of  a  segment  of  poor  conductivity  more 
closely  matches  the  observed  |E‘  x|  variation  than  the  elevation  effect.  The  models 
demonstrate  that  irregularities  in  both  profile  and  earth  conductivity  play  significant 
roles  in  determining  the  propagation  of  Mf'  over  ridges. 


10.  Berry,  L.  A.  (1967)  Radio  propagation  over  a  Gaussian-shaped  ridge,  IEEE 

Trans.  Antennas  Propag.  ^5(No.  5):701-7  02. 

11.  FAO-Unesco,  Soil  Map  of  the  World,  Vol.  II,  Unesco-Paris,  1974. 
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GAUSSIAN  RIDGE  AT  ifSa 
(h  •O.Akm,  w  •  T.Skm) 


Ki^ure  10.  'I  he  AUcnuation 
l''uncti<)n  ,\V  |  )n  Iho  Surfaci' 
)f  ilii'  Hidge  Shown  al  the 
Bottom  ol’  the  I'igure,  lor 
I'wo  \  alues  of  l^arth 
Conduct  ivity 
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Figure  11.  The  Attenuation 
Function  fwl  jn  the  Surface 
of  the  Inhomogeneous  Path 
Shown  at  the  Bottom  of  the 
Figure,  for  Three  Values 
of  a 
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6.  I'K.NKTRATIO.N  IIF.LOW  TIIK  l)FSI.RT  FI.OOR 


Attonuation  of  .Ml'  t'ii-lds  with  depth  in  homogeneous,  well-eondueting  soil  is 
described  bv 


-v/6  -jv/6 

H  ^  II  e  '  e  ‘ 
o 

wlu'i'e  is  the  magnitude  of  the  field  at  the  surface,  y  is  the  depth  below  the  sur¬ 
face  and  6  ,  is  the  depth  of  penetration  or  skin  depth.  Attenuation  measurements  of 
tne  magnetic  fields  penetrating  the  desert  soil  were  made  by  lowering  a  loop  an¬ 
tenna  into  18-in.  diameter  holes  drilled  approximately  6-ft  deep  by  a  commercial 
power  auger,  and  recording  the  decibels  chang('  from  the  surface  reading.  An 
effective  skin  depth  was  then  calculated  via 

°s  [  dl3]  -  III  dBj  • 

A  total  of  30  holes  were  drilled  at  various  locations  as  indicated  in  l-'igure  I.  The 
relative  attenuation  was  re;id  to  within  tibout  ±  0.  25  dB  on  the  .N1V1-25T  panel  meter, 
and  the  depth  was  measured  accurately.  The  average  skin  depth  calculated  was 
28.  25  ft  with  a  standard  deviation  of  13.  00  ft. 


7.  DISCb.S.SIOiN  AiNU  COXlLUSIOiN.S 

The  data  in  this  report  can  be  used  to  estimate  total  propagation  path  loss  for 
a  point-to-point  communication  link  between  terminals  buried  in  the  desert  south¬ 
east  of  Yuma  and  separated  by  mountain  ridges.  The  desert  valley  soils  have  rela¬ 
tively  high  apparent  conductivity,  with  a  correspondingly  low  attenuation  of  propa¬ 
gating  MF  radio  waves.  The  mountain  ridges  can  have  a  large  effect  which  depends 
on  the  location  of  the  terminal.  The  mountains  can  be  modeled  to  predict  the 
attenuation,  but  good  accuracy  requires  that  both  the  terrain  profile  and  the  ground 
constants  be  known.  An  overburden  of  desert  soil  will  result  in  additional  loss. 

The  measurements  and  portions  of  the  modeling  work  in  this  report  were  pre¬ 
sented  to  an  MX  C  working  group  at  a  technical  interchange  meeting  at  Norton  AFB, 
California  on  14  December  1977. 
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